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Based <>n purely matrix-algebraic arguments, we prove three new results on the eigenvalue 

problem [Zimm, L956] arising from modeling the hulk hydrodynamic and dielectric properties of very 
dilute polymer solutions: 

1. Let V he the number of identical segments joining /V+l heads of a head/spring model for any 
polymer molecule. Then for any N, the well-known hydrodynamic interaction matrix H is 

positive definite if the so-called interaction parameter h* is less than (2v2 + l)/7 or 0.547. 

2. Let A he the tridiagonal matrix governing the elastic link force on each head. Then il H is 

positive definite, there exists a matrix (^ such that (J -1 H Ay = (j 7 AQ = A, where A is a 

diagonal matrix. 
'.]. (Equivalence Statement due to Lodge and Wu [1971]). The nonzero eigenvalues of the matrix 
product H A are completely determined hy solving the eigenvalue problem associated with a 

symmetric matrix S, where Sy = H u + Hi- it j-i — fff-u - Hu i» l »7 = 1*2 N. 

To illustrate the significance of these results, numerical output for N as large as 300 based on a 
Fortran program for several values of h* is given. 

Key words: Bead/spring model; dilute polymer solutions; eigenvalue; hydrodynamic interaction: 
matrix algehra; necklace model; polymer physics. 



1 . Statement of Problem 

Some twenty years ago, Zimm [l] 1 formulated a linear, second-order partial differential equation 
for a distribution function \fj which depends on time and 3(N+1) coordinates xq, y& Zq, . . . , x N , y N , 
z N , of the N + 1 beads, for modeling the bulk behavior of very dilute polymer solutions under the 
influence of external force, Brownian motion, and hydrodynamic interaction among the beads of the 
necklace model. The mechanical model for each polymer molecule is that of a chain of N identical, 
ideally elastic segments joining /V+l identical beads with complete flexibility at each bead. Two 
length parameters are of interest in this model: a h , the so-called hydrodynamic radius of the bead, 
and b , the root mean square of the segment length. The ratio a of the two length parameters (a = 
cifjbo), and the number N of elastic segments completely characterize the mathematical problem 



*On leave during 1975-76 as a ComSci Fellow, Nuclear Regulatory Comm., Wash., D.C. 2055! 
1 Figure in brackets indicate the literature references on page 281. 
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through the following two symmetric, square matrices of order 7V+1: 



(1.1) 
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(1.2) 



H 



1 a a.2"* a.3~* a.N~l 

a 1 a a.2-t a(W-l)"* 

a. 2~* a 1 a a(N-2)^ 

a./V-* a(/V-ir* 1 



The symmetry of the two matrices, A and H , and an assumption based on physical arguments that 
the eigenvalues of the matrix product H A are distinct, led to an eigenvalue problem with a 
coordinate transformation matrix y such that a similarity transformation with (J diagonalizes H_A, 
and congruent transformations with (J diagonalize A_and H_ as follows: 



(1.3) 
(1.4) 

(1.5) 



Q-' HA Q 


= a , 




y T t y = 


= M , 




y-'H [q-t = 


= N = 


= A M" 1 



where A = diag(A , Xi, . . ., \\), M = diag(/x , /x x , . . ., /jl n ), and N = diag(^ () , v u . . ., v N ), are 

diagonal matrices, and the superscripts _1 and t in eqs (1.3) through (1.5) denote the inverse and 
the transpose of a matrix, respectively. A closed form solution for a specific linear laminar flow with 
a transverse gradient constant over all beads, was first given by Peterlin and appeared at the end of 
Zimm's paper [1, p. 278] as a note added in proof. Computationally speaking, the modeling problem 
is reduced to the construction of the transformation matrix (J such that both A and M_ are known 

for constructing the distribution function i//. 

The purpose of this paper is to derive some new results based on matrix algebraic arguments 
and another obvious property of the matrix H, namely, its positive-definiteness under certain 

conditions, which Zimm [1] failed to exploit. In section 2, we first show that H is positive definite if 

the parameter a is positive and is bounded from above by a constant a ma x which equals 1.0 for N = 
1, 0.872 for W = 2, 0.815 for W = 3, and 0.773 for /V > 4. In section 3, we prove our main result 
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that for positive definite H, there exists some Q such that A = M and N = J_, i.e., kj = /Xj and Vj 

= 1 for j — 0, 1, 2, . . ., N. In section 4, we furnish an algebraic proof as an alternate to the proof 
of a result due to Lodge and Wu [2], which we subsequently use in calculating the eigenvalues of 
the matrix product H A. In section 5, we discuss the significance of our three mathematical results. 

To illustrate the significance, we include in an appendix a sample Fortran program and output for 
several values of the nondimensional parameter a and for N as large as 300. 



2. Positive Definiteness of the Interaction Matrix 

As shown in (1.2), the so-called interaction matrix H is not only symmetric, but also a "stripe 

matrix, " i.e., all entries along a line parallel to the main diagonal are equal to one another. 
Furthermore, H depends on a scalar parameter a, which is required to be positive (a > 0) from 

physical considerations. 2 We now wish to show I hat, for any given order of the matrix H, there is an 

upper bound a max for the parameter a such that H is positive definite for < a < a max . 

Since H is a stripe matrix, an application of a standard result in matrix algebra 3 Implies the 

corollary that H is positive definite, if and only if, the determinant of H, i.e., (let II_, is positive for a 

specific order IS + 1 as well as for all orders less than /V 4- 1. It is convenient to introduce a new 
matrix B and a new positive constant c as defined below: 



(2.1) 



(2.2) 



(i 



B = 



i 

c 1 2"i 3"* W-4 

\ c 12"* .. (ZV-l)-i 

2"* 1 c 1 (/V-2) * 

N-i (/V-ir* (/V-2)-* (/V-3)-* c 



Since det B = c N+l det H, c>0, our problem is reduced to the study of the positivity of det B. 

From the property of a determinant that its value is unaltered if any multiple of one line is 
added to a different, parallel line, 4 let us perform the following operations on the matrix B: 

Step I Keep the first row of B fixed, and replace the ;th row with the difference between they'th 
and the (7 — Dth rows for; = 2,3, . . ., N + l. 

Step 2 Keep the first column of the new matrix following step 1 fixed, and replace the jlh column 
of the new matrix with the difference between the jth and the (/ — l)th columns of the new 
matrix for) = 2,3, . . ., W + l. 



2 We have already seen that hy definition, a{=a,Jb {> ) is required to he positive ii the hydrodynamic radius a,, is positive. In Zimm's formulation [1, pp. 270— 
271], the parameter is defined hy the expression (ott'i * p tj '6 ', which is positive for p > 0, 17 > 0, and b > 0. 
" see, e.g., riotin [4, pp. fttf-ZOUJ. 
4 See, e.g., Hohn [3, p. 37]. 
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The resulting matrix, to be denoted by B', is again symmetric as shown below: 



(2.3) 



B' = 



c 

1-c 

2-4-1 

3"*-2-* 



1-c 
2(c-l) 

-2-H2-c 
-3"* + 2(2)-*-l 



N-i-{N-l)-i -A~*+2(A - i)-i-(7V-2)-i 



2-i-l . 
-2-H2-C. 

2(c-l) . 
-2"*+2-c. 



Since det B' = det B, we have further reduced the study of the positive definiteness of the 

interaction matrix H to a study of the matrix B_\ As it will become clear very shortly, our 

derivation so far is intended to pave the way for the application of a sufficient condition for the 
positive definiteness of H or B'. 

We now observe that for the lowest order of H, B, or B\ i.e., all the 2x2 matrices given by 

(1.2), (2.2), or (2.3), the positivity of det H (and thus det B^or det ET) always yields the condition 

that either a < 1 or c > 1. It follows that the symmetric matrix B' has the property that all its 

diagonal elements are positive. A well-known theorem in matrix algebra states that a symmetric 
matrix with positive diagonal elements is positive definite if the matrix is "strongly diagonal 
dominant," i.e., the absolute value of the diagonal element at every row exceeds the sum of the 
absolute values of all off-diagonal elements in that row: 5 



(2.4) 



\B'„ 



N 

> 2 |fl*l 

fc=o,i,.. 

(k=j) 



for; = 0,1,2, . . ., N. 



We denote the generic of B' as B jk ,j,k = 0,1,2, . . ., N. 

To apply (2.4), we consider two ranges of the values of c for which one can write down easily 
the A^ -hi inequalities: 



Case I (c > 2 - 2"*): 

(2.5) ;=0: c > c - N~*. 

(2.6) ;=1: 2(c - 1) > 2(c - 1) - (A-l)-± + N~K 

(2.7) 7=2: 2(c - 1) > 2(c - 1) - (A-2)-* + (IV- 1)^. 

(2.8) j= A: 2(c - 1) >c - 1. 
Case 2 (1 < c < 2 - 2"*): 

(2.5) 7=0: c >c - N~K 

5 See. e.g.. Schwaiz [4, pp. 17-19]. 
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(2.9) ;=1: 2(c - 1) > 1 - 2"* 

(2.10) or, 2(c - 1) > 2 - 3(2)-* + 3"* 

(2.11) or, 2(c - 1) > 2 - 2(2)-* - (/V-l)~* + /V* 

(2.12) ; = 2: 2(c - 1) > 3 - 2(2)"* - c 

(2.13) or, 2(c - 1) > 5 - 3(2)-* - 2c 

(2.14) or, 2(c - 1) > 6 - 5(2)-* - 2c + 3"* 

(2.15) or, 2(c - 1) > 6 - 4(2)"* - 2c - (/V-2)"* + (N-l)"* 

(2.12) ;=3: 2(c - 1) > 3 - 2(2)-* - c 

(2.13) or, 2(c - 1) > 5 - 3(2)"* - 2c 

(2.14) or, 2(c - 1) > 6 - .5(2)-* - 2c + 3"* 

(2.15) or, 2(c - 1) > 6 - 4(2)~* - 2c - (/V-3)"* + (/V-2)-* 



(AT = 2), 

(# = 3), 

(/V i= 4). 

(tf = 2), 

(N = 3), 

(JV = 4), 

(/V ^ 5). 

(JV = 3), 

(iV - 4), 

(N = 5, 

(/V i= 6). 



(2.16) y = /V: 2(c - 1) > 3 - 2(2)"* - c 



(N ^ 2). 



What distinguishes the above two cases is the fact that for case 1, every off-diagonal element of B' 

is negative, whereas for case 2, care must be exercised on elements B' 2 & B' 3 % B' 34 . . ., B' NtN _ u each 
of which equals — 2~* + 2 — c. An examination of inequalities (2.5) through (2.8) for case 1 implies 
that they are always satisfied for every N. On the other hand, for case 2, the controlling inequalities 
are given by (2.9), (2.10), and (2.11). Combining all results derived in this section, including (2.1) 
where a = c _I , we obtain various upper bounds a max f° r tne parameter a for various /V as listed in 
table 1. 



Table 1. Upper bounds for a or h* for positive definite H 



Number of Chains 


« max 




"max — * u*niax' 


(Exact) 


(Approx.) 


(Approx.) 


N m 1 


1 




1.0 


0.707 


N = 2 


^(12 + 2\/2) 




0.872 


0.617 


N = 3 


3216 + 1278V2 - 484V3 - 
4273 


- 288V6 


0.815 


0.576 


N = 4 and above 


y(4 + V2) 




0.773 


0.547 
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We have included in table 1 a list of upper bounds for a parameter h* which appears frequently in 
the literature with h* = 2~* a. Since most of the applications of the bead/spring model require more 
than four chains (JV > 4), the key result, < h* < 0.547, a criterion which assures positive 
definiteness of H, will be used in the next section to achieve a significant simplification of Zimm's 

eigenvalue problem. 



3. Diagonalization — A Simplification of Zimm's Result 

The derivation of an upper bound for the parameter a (0 < a < 0.773) or, equivalently, the 
parameter A* (0 < h* < 0.547), with h* = 2~* a, as shown in the last section, enables us to inject 
the positive definiteness of the matrix H as an additional matrix property into Zimm's eigenvalue 

problem. We begin with another standard result in matrix algebra that any symmetric, positive 
definite matrix can be represented by the square of another symmetric matrix. 6 Hence, for < A* 
< 0.547, the symmetric matrix H is positive definite and there exists a symmetric C such that H = 

C 2 , C T = C. Since the matrix A as defined in (1.1) is symmetric, the matrix product C A C is 

also symmetric. By the well-known diagonalization theorem on symmetric matrices, 7 we assert that 
there exist an orthogonal U and a diagonal D^ such that 

(3.1) C A C = U" 1 D U; U 7 U = 1. 

By choosing Q as the matrix product C U "V eqs (1.3) through (1.5) can be simplified as follows: 

(3.2) Q- 1 H A y = U C^C 2 ^ CU 1 =D =A, 

(3.3) Q T A (J = U C A C U 1 = D = M 

(3.4) Q- 1 H(Q- 1 ) T = UC-'^C-'U- 1 = 1 = N. 

Thus we have shown that for < h* < 0.547, there exists a transformation matrix 
Q such that A = M, and N = jL 

4. Equivalence of Two Eigenvalue Problems — An Algebraic Proof 

For the bead and spring model of dilute polymer solutions, the result of last section shows that 
it is sufficient to find the /V+l eigenvalues of the matrix product H A if the criterion < a < 0.773 

or < h* < 0.547 is satisfied. Since A is singular and all eigenvalues are distinct, det (H A) 

vanishes and exactly one of those /V + 1 eigenvalues must be zero, i.e., \ — 0. 

To find the other A^ eigenvalues, \ u . . ., k N , Lodge and Wu [2] introduced an equivalence 
statement which asserts that the eigenvalue problem for the unsymmetric product H A of order /V + 

1, aside from the addition of a zero eigenvalue, is equivalent to that of a symmetric matrix 



"See. e.g., Halmos [2. p. 139]. 

7 See, e.g.. Halmos [2. p. 157. bottom lines]. 
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S of order N, where the elements S t j of S are combinations of elements of H as follows: 



(4.1) 



Si j = H (j + //,-i.j-i - //«-i,j - ff|.j-i, 



i,j=lA 



N. 



The equivalence statement is important because the matrix S for the reduced problem is symmetric 

and, therefore, computationally more attractive than the matrix product H A. To prove that 

statement, Lodge and Wu [2] introduced a nonsingular matrix A € with elements (A e )ij = A ti for 

each ij = 0,1, . . ., N, except for i = j = where (A € ) 00 — Aqq + e 2 . Using a theorem 8 on the 
continuity of eigenvalues, a triangular decomposition of the matrix A € , and a property of the 

determinant of a matrix product, it was shown that as e ^ 0, the nonzero eigenvalues of H A e 

approach those of H A and S. 

The purpose of this section is to show that Lodge and Wu's equivalence statement is true by a 
purely algebraic argument without resorting to the continuity theorem involving the parameter e. We 
first observe that there exists a triangular decomposition of A such that A = L 7 L: 



(4.2) 



L = 












1 


1 








-1 


1 



Let A, be any nonzero eigenvalue of H A, i = 1,2, . . ., N. Let x, be the corresponding eigenvector 
of Xj. By definition, we have 



(4.3) 



H A x j = kjXj. 



Substituting A = L 7 L into (4.3), we obtain 



(4.4) 



H L 7 L x t = X/Xj. 



By definition, x t is a nonzero vector satisfying (4.3), hence the right-hand side of (4.4) cannot vanish. 
This implies that L x t is a nonzero vector (proof by contradiction). Therefore it makes sense to 

premultiply both sides of (4.4) by the singular matrix L: 



(4.5) 



L H L 7 Lx t = ki L %i. 



Denoting L *, by a nonzero vector y h we obtain 



(4.6) 



L H L 7 >,= \,y„ 



* Lodge and Wu [2] attribute this theorem to Ostrowski [6]. 
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Equations (4.3) and (4.6) show that the nonzero eigenvalues of H A and L H L r are identical. Some 
algebraic manipulation with L HL T , which is symmetric, shows that its nonzero eigenvalues are 
given by those of S with components S u defined in (4.1). 



5. Significance of Results 

To discuss the significance of our mathematical results, it is appropriate here to introduce a 
closed-form solution of the distribution function 0, originally due to Peterlin [1, p. 278 footnote], for 
a specific linear laminar flow with a constant transverse velocity gradient k expressible through a 
nondimensional gradient parameter /3: 9 



(6.1) 
where 

(6.2) 



*/> = */>! </>2 $i $N* 



ipi = C i exp 



Mi 



2 1 4- £ 2 /\? 



& 






2/3* 



+ d + -fr)£? 



(l + 



tf 



vt 



C i = (jjii/ff) 312 (1 + /3 2 /k*) 1/2 , and \ 4 , /x 4 , i = 1,2, . . ., N, are the nonzero diagonal elements of 
matrices A and M as defined in (1.3) and (1.4), respectively. The symbols g u £ t -, 77 ,- in (6.2) denote 

the nondimensionalized transformed coordinates of the (i + l)th bead as defined below: 



(6.3) 



(6.4) 



j=i 

v 



(6.5) 



T?i= I 6c 1 [(>-']«*). 



The principal result of this paper has been the demonstration that the diagonal elements of both 
matrices A and M are identical for each i, i — 0,1,2, . . ., N, if the hydrodynamic interaction 

parameter h* is less than 0.547. This means that once we have computed the values of A,, i = 0,1,2, 
. . ., N, from the eigenvalue problem of the matrix product H A as defined in (1.3), it is no longer 

necessary, as understood previously by Zimm [1], to compute Q explicitly in order to calculate the 
diagonal elements fi h i = 0,1,2, . . ., N, of the matrix M as defined in (1.4). 

The result A, = fi h i = 0,1,2, . . ., A, for fe* < 0.547 is also significant in discrediting some 
computational results in the literature such as Zimm, Roe, and Epstein [7], etc., where the 
transformation matrix (J, being not explicitly known, was expeditiously approximated by the 

orthogonal transformation matrix Qr of the well-known Rouse model of zero hydrodynamic 

interaction (i.e., H = I). Under that circumstance, the matrix M is given by Q^i A Q_ r or 



9 The flow field is specified by the velocity vector v = Ik y. 0, 0). and the parameter fi is related to k by the' formula f3 
Boltzman's constant and T is the absolute temperature. 

10 See. e.g.. Zimm, Koe and Epstein [7], Hearst [8], and Tschoegl [9]. 
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Q.R AQfo an< ^ ^ e eigenvalues /jlj are given by 4 sin 2 \j tt/2{N+1)] which contradict our exact result 

with \j = n h and M _1 A = N=]_as shown in (3.4). 

Computationally speaking, we are interested in the solution of the eigenvalue problem (4.3), for 
large /V, and particularly, in the smaller values of A*. The equivalence statement of Lodge and Wu 
[2] provides the basis for computing those eigenvalues very economically with the aid of a large-scale 
computer. The exact computer-generated calculation renders many earlier approximation schemes 
obsolete. 10 Without the benefit of the equivalence statement, Thurston and Morrison [10] used a 
computer to obtain the eigenvalues of the matrix product H A for values of N up to 15. Based on 

their equivalence statement, Lodge and Wu [11] were able to extend the computational capability 
and reported some numerical results for bead/spring models with as many as 300 links. 
Unfortunately, Lodge and Wu did not exhibit all eigenvalues, nor the computer program with which 
eigenvalues could be generated for interesting applications. In the appendix we list a Fortran 
program and output to illustrate the significance of this paper and to enable others to extend these 
calculations, if they so desire. 



We are indebted to Morris Newman, Senior Research Mathematician, National Bureau of 
Standards, Washington, D.C., for substantial original contribution and critical comments during the 
course of this investigation. We also wish to thank Martin Cordes of NBS Applied Mathematics 
Division for technical consultation related to the computational aspect of this work. 
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7. Appendix 



The following is a Fortran V driver program as executed on a UNIVAC-1 108 Exec 8 computer 
at National Bureau of Standards Laboratory at Gaithersburg, Md. To simplify our programming 
effort, we take advantage of the well-documented subroutines furnished by the International 
Mathematical and Statistical Libraries, Inc., of Houston, Texas. For a complete description of these 
subroutines used in our program, the reader may consult the reference manual entitled 'TMSL 
Library 2, Edition 4 (FORTRAN V), UNI VAC 1100 Series, Reference Manual, November 1974." 
The appearance of trade names such as UNIVAC and IMSL does not imply an endorsement of their 
product or services. 
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1* 

2* 




COMMON A(300),B(300),D(300),E(300),E2(300),HMINOR(45150) 

N = 300 


3* 




HSTAR = 0.4 


4* 




DO 20 1 = 1,N 


5* 




A(I) = HSTAR* SQRT(2./FLOAT(I)) 


6* 


20 


CONTINUE 


7* 




B(I) = 1. - A(l) 


8* 




DO 30 1 = 2,N 


9* 




B(I) = A(I - 1) - A(I) 


10* 


30 


CONTINUE 


11* 




DO 50 I = 1,N 


12* 




II = I * (I+D/2 


13* 




HMINOR(II) = 2. * B(I) 


14* 


50 


CONTINUE 


15* 




Nl = N - 1 


16* 




DO 60 1 = 1, Nl 


17* 




N2 =.N - I 


18* 




DO 70 J = 1, N2 


19* 




K = ((I+J)*(I+J-l)/2) + J 


20* 




HMINOR(K) = B(I+1) - B(I) 


21* 


70 


CONTINUE 


22* 


60 


CONTINUE 


23* 




CALL EHOUSS (HMINOR, N, D. E, E2) 


24* 




CALL EQRT2S (D, E, N, DUMMY, 1, IERR) 


25* 




IFdERR.NE. 0)GOTO 100 


26* 




SUMD1 = 


27* 




SUMD2 = 


28* 




DO 80 I = 1,N 


29* 




SUMD1 = SUMD1 + l./ENI) 


30* 




SUMD2 = SUMD2 + l./D(I)/D(I) 


31* 


80 


CONTINUE 


32* 


200 


FORMAT (1H1, 10X, 4HN = , 14. 20X 8HHSTAR = ,F6.2,//) 


33* 




WRITE (6,200) N, HSTAR 


34* 


250 


FORMAT (1H, 4\. 42HNONZERO EIGENVALUES FOR MATRIX PRODUCT H*A,//) 


35* 




WRITE (6.250) 


36* 


320 


FORMAT (1H0, 10F12.5) 


37* 




WRITE (6,320) (D(I). 1=1, N) 


38* 


330 


FORMAT (1H0,8X,E12.5,2X,20H(SUM OF RECIPROCALS)) 


39* 




WRITE (6,330) SUMD1 


40* 


340 


FORMAT (1H0,8\.E12.5.2X.27H(SUM OF RECIPROCAL SOUARES)) 


41* 




W RITE (6,340) SUMD2 


42* 


100 


W RITE (6,400) IERR 


43* 


400 


FORMAT (1H1. 10X. 36HROUTINE FAILS AT EIGENVALUE NUMBER .14) 


44* 




END 
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For N = 300, and h* = 0.4 (corresponding to a = 0.566), the following output was obtained at a 
total run cost of $47.91 for a total run time of 5 min 42 s. The CPU cost was $24.21 with a run time 
of 5 min 36 s. 



N = 300 HSTAR = .40 

NONZERO EIGENVALUES FOR MATRIX PRODUCT H*A 



.11539-02 
.45574-01 
.11251 + 00 
.18976 + 00 
.27226 + 00 
.35690+00 
.44150+00 
.52447 + 00 
.60458+00 
.68091 + 00 
.75277+00 
.81967+00 
.88130+00 
.93748+00 
.98816+00 
.10334+01 
.10734+01 
.11083 + 01 
.11385 + 01 
. 11643+01 
.11860+01 
.12041+01 
.12190+01 
.12310 + 01 
.12435 + 01 
.12480+01 
.12536 + 01 
.12577 + 01 
.12600 + 01 
.12670+01 



.35945-02 
.51576-01 
.11988+00 
.19784+00 
.28067+00 
.36539+00 
.44990+00 
.53263+00 
.61240+00 
.68830+00 
.75969+00 
.82608+00 
.88716+00 
.94279+00 
.99293+00 
.10376+01 
.10771 + 01 
.11 116+01 
.11413+01 
.11666+01 
.11880+01 
.12058+01 
.12203+01 
.12321 + 01 
.12414+01 
.12486+01 
.12541 + 01 
.12580+01 
.12606+01 
.12621 + 01 



.67180-02 
.57762-01 
.12734+00 
.20596+00 
.28909+00 
.37388+00 
.45827+00 
.54075+00 
.62017+00 
.69566+00 
.76656+00 
.83243+00 
.89297+00 
.94806+00 
.99765+00 
.10418+01 
.10808+01 
.11147 + 01 
.11440+01 
.11689+01 
.11899 + 01 
.12073 + 01 
.12216+01 
.12331 + 01 
.12422 + 01 
.12492+01 
.12545+01 
.12583+01 
.12608+01 
.12622+01 



.10383-0] 

.64119-01 
.13489+00 
.21412+00 
.29754+00 
.38236+00 
.46663+00 
.54885+00 
.62791 + 00 
.70296+00 
.77338+00 
.83872+00 
.89873+00 
.95326+00 
.10023 + 01 
.10460+01 
.10844+01 
.11179+01 
. 1 1467 + 01 
.11712 + 01 
.11918+01 
.12089+01 
.12229+01 
.12341 + 01 
.12430+0] 
.12492+0] 
.12550+01 
.12586 + 01 
.12610+01 
.12623+01 



.14480-01 
.70634-01 
.14252+00 
.22233 + 00 
.30599+00 
.39084+00 
.47496+00 
.55691 + 00 
.63531 + 00 
.71022+00 
.78015+00 
.84497+00 
.90443+00 
.95841 + 00 
.10069+01 
.10500+01 
.10880+01 
.11210+01 
.1 1493 + 01 
.11734+01 
.11937+01 
.12104+01 
.12241 + 01 
.12351 + 01 
.12438+01 
.12504+01 
.12554+01 
.12589 + 01 
.12612 + 01 
.12674+01 



.18953-01 
.77297-01 
.15023+00 
.23058+00 
.31446+00 
.39932+00 
.48328+00 
.56495+00 
.64326+00 
.71743+00 
.78686+00 
.85116+00 
.91008+00 
.96351 + 00 
.10115+01 
.10541 + 01 
.10915 + 01 
.1 1210 + 01 
.11519 + 01 
.11756+01 
.11955+01 
.12119+01 
.12253 + 01 
.12361+01 
.12445 + 01 
.12510+01 
.12558+01 
.12592 + 01 
.12614+01 
.12624+01 



.23749-01 
.84098-01 
.15801+00 
.23886+00 
.32294+00 
.40778+00 
.49157+00 
.57294+00 
.65088+00 
.72459 + 00 
.79353+00 
.85730+00 
.91567+00 
.96855+00 
.10160+01 
.10580 + 01 
.10949+01 
.11270 + 01 
.11545+01 
.11778+01 
.11973 + 01 
.12134+01 
.12265 + 01 
.12370+01 
.12452 + 01 
.12515+01 
.12562+01 
. 12595 + 01 
.12615 + 01 
.12625 + 01 



.28837-01 
.91028-01 
.16586+00 
.24717+00 
.33142 + 00 
.41623+00 
.49983+00 
.58091 + 00 
.65845 t 00 
.73171 + 00 
.80014+00 
.86338 + 00 
.92120+00 
.97353+00 
.10204+01 
.10619+01 
.10984+01 
.11299+01 
.11570+01 
.11799 + 01 
.11991 + 01 
.12148+01 
.12277+01 
.12379 + 01 
.12460+0] 
.12521 + 01 
.12566+01 
.12597+01 
.12617 + 01 
.12625 + 01 



.34185-01 
.98078-01 
.17377+00 
.25551+00 
.33991+00 
.42467+00 
.50807+00 
.58883 + 00 
.66598+00 
.73878 + 00 
.80670 + 00 
.86941+00 
.92668+00 
.97847 + 00 
.10248+01 
.10658+01 
.11017+01 
.11328+01 
.1 1595 + 01 
.1 1820 + 01 
.12008+01 
.12163 + 01 
.12288+01 
.12388+01 
.12466 + 01 
.12526 + 01 
.12570 + 01 
.12600 + 01 
.12618+01 
.12625 + 01 



.39772-01 
.10524+00 
.18174+00 
.26387+00 
.34840+00 
. 13309+00 
.51629 + 00 
.59673+00 
.67347+00 
.74580+00 
.81321+00 
.87538+00 
.93211+00 
.98334+00 
.10291+01 
.10696+01 
.11051+01 
.11357 + 01 
.11619 + 01 
.11840+01 
.12025+01 
.12176 + 01 
.12299 + 01 
.12397+01 
.12473+01 
.12531+0] 
.12573+0] 
.12602 + 01 
.12619+01 
.12625+01 



By changing the first three statements of the original program, one can find the eigenvalues for any 
TV and A* as long as the core storage is not exceeded and the value of /?* is less than 0..547. To 
illustrate this point, we show below for N = 250 and /z* = 0.3 the first three statements of the 
modified program and the appropriate output. The total cost for the run was $25.29 with a total run 
time of 3 min 19 s. The CPU cost was $13.96 and the CPU time was 3 min 14 s. 



COMMON A(250),B(250),D(250),E(250),E2(250),HMINOR(31375) 

N = 250 
HSTAR = 0.3 
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N = 250 HSTAR = .30 

NONZERO EIGENVALUES FOR MATRIX PRODUCT H*A 



.11737-02 
.48849-01 
.12499 + 00 
.21726+00 
.32025 + 00 
.43042+00 
.54511 + 00 
.66213 + 00 
.77969 + 00 
.89622 + 00 
.10104+01 
.11211+01 
.12272+01 
.13280+01 
.14228+01 
.15108+01 
.15916+01 
.16648+01 
.17301+01 
.17871 + 01 
.18356+01 
.18755+01 
.19067+01 
.19290+01 
.19424+01 



.36804-02 
.55504-01 
.13360+00 
.22714+00 
.33099+00 
.44173+00 
.55674+00 
.67389 + 00 
.79141 + 00 
.90776+00 
.10216+01 
.11319+01 
.12376+01 
.13378+01 
.14319+01 
.15192+01 
.15993+01 
.16717+01 
.17362+01 
.17923+01 
.18400+01 
.18790+01 
.19093+01 
.19308+01 
.19433+01 



.69245-02 
.62402-01 
.14237+00 
.23713+00 
.34180+00 
.45308+00 
.56039+00 
.68565+00 
.80312 + 00 
.91928+00 
.10328+01 
.11427+01 
.12478+01 
.13475+01 
.14409+01 
.15276+01 
.16069+01 
.16785+01 
.17421+01 
.17975+01 
.18443+01 
.18825+01 
.19119+01 
.19324+01 
.19440+01 



.10767-01 
.69531-01 
.15128+00 
.24721 + 00 
.35268+00 
.46446+00 
.58006+00 
.69742+00 
.81482 + 00 
.93077+00 
.10440+01 
.11534+01 
.12581 + 01 
.13571 + 01 
.14499+01 
.15358+01 
.16144+01 
.16853+01 
.17481 + 01 
.18025+01 
.18485+01 
.18858+01 
.19143+01 
.19340+01 
.19447+01 



.15099-01 
.76877-01 
.16033 + 00 
.25739+00 
.36362+00 
.47589+00 
.59174+00 
.70918+00 
.82650+00 
.94224+00 
.10551 + 01 
.11641+01 
.12682+01 
.13667+01 
.14588+01 
.15440+01 
.16218+01 
.16919+01 
.17539+01 
.18075+01 
.18526+01 
.18891 + 01 
.19167+01 
.19355 + 01 
.19453+01 



.19865-01 
.84432-01 
.16951 + 00 
.26766+00 
.37462+00 
.48735+00 
.60345 + 00 
.72094+00 
.83817+00 
.95368+00 
.10662 + 01 
.11748+01 
.12783+01 
.13762+01 
.14676+01 
.15521 + 01 
.16292+01 
.16905+01 
.17596+01 
.18124+01 
.18567 + 01 
.18922+01 
.19190+01 
.19368+01 
.19458+01 



.25015-01 
.92184-01 
.17882 + 00 
.27802+00 
.38568+00 
.49894+00 
.61516+00 
.73270 + 00 
.84982 + 00 
.96508+00 
.10773+01 
.11854+01 
.12884+01 
.13856+01 
.14764+01 
.15602+01 
.16365+01 
.17050+01 
.17653+01 
.18172+01 
.18606+01 
.18953+01 
.19212+01 
.19381+01 
.19462 + 01 



.30516-01 
.10012 + 00 
.18826+00 
.28846+00 
.39679 + 00 
.51037+00 
.62689 + 00 
.74446+00 
.86145+00 
.97646+00 
.10883+01 
.11959+01 
.12984+01 
.13950+01 
.14851 + 01 
.15682+01 
.16437+01 
.17114+01 
.17709+01 
.18220+01 
.18645 + 01 
.18983+01 
.19233+01 
.19393+01 
.19465+01 



36336-01 
,10824+00 
19781+00 
29898+00 
.40795+00 
.52192+00 
,63863+00 
,75621+00 
87306+00 
,98780+00 
.10993+01 
,12064+01 
,13083+01 
,14043+01 
.14938+01 
.15761 + 01 
,16508+01 
,17177+01 
.17764+01 
.18266+01 
.18683+01 
.19012+01 
,19253+01 
.19405+01 
,19467+01 



.42455-01 
.11653+00 
.20748+00 
.30958+00 
.41916+00 
.53350+00 
.65038+00 
.76795+00 
.88465 + 00 
.99912 + 00 
.11102+01 
.12168+01 
.13182 + 01 
.14136+01 
.15023+01 
.15839 + 01 
.16579+01 
.17239+01 
.17818+01 
.18312+01 
.18719+01 
.19040+01 
.19272+01 
.19415 + 01 
.19469 + 01 



(Paper 80A2-443) 
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